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Decay amplitudes of D(D S ) — > /o(980)X, X = tt, K, are compared to experimental 
branching ratios with the aim of singling out the poorly known D — > / (980) transition 
form factor in these amplitudes. Since the other elements of the amplitudes are either 
calculable in an effective QCD theory using operator product expansion or are known 
from experiment (e.g. the pion and kaon decay constants), we can take advantage of these 
reactions to constrain the transition form factors obtained in relativistic quark models £Q . 
In these models, the /o(980) wavefunction requires an unknown size parameter for both 
its non-strange uu(dd) and strange ss components, which we fit to the D(D S ) decay data. 

1. MOTIVATION 

Scalar mesons have been a recurrent topic over the past 30-40 years and while the a has 
been a longstanding open question since the 1960s, the / (980) and a (980) were firmly 
established in tttt scattering experiments in the 1970s 0. The constituent structure of 
the scalar mesons, nevertheless, has been and still is a controversial issue. For a general 
overview we refer, for instance, to the Particle Data Group review [3] and references 
therein. The known ++ mesons fall into two classes: near and about 1 GeV and in the 
region 1.3 — 1.5 GeV. The scalar objects below 1 GeV form an SU(3) flavor nonet. This 
nonet contains two isosinglets, an isotriplet and two strange isodoublets. Among these 
lighter scalars, the isosinglet /o(980) and the isotriplet a(980) are fairly narrow: T = 
40 — 100 MeV. Both couple strongly to the KK channel and lie close to the KK threshold 
at 987 MeV. This closeness to threshold distorts the shape of their resonant structure and 
the description of the f and a requires a coupled-channel scattering analysis. 

The simple quark model views these scalar mesons as orbitally (L = 1) excited qq states 
and has been advocated by Tornqvist and Roos 4J. However, some studies jS] tend to 
favor four-quark configurations of the scalar mesons and so do coupled-channel analyses [H] 
or potential models of molecular states strongly coupled to tttt and KK thresholds [7j. 
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In any case, the precise four-quark picture is still open to debate since a coupled-channel 
model implies two qq pairs while the configuration q 2 q 2 requires two di-quarks. 

In the past few years, scalar mesons have emerged in three-body decays of heavy mesons 
such as D — > 7T7T7r, ttttK, KKK measured at CESR, DESY and Fermilab [8-16] and B — > 
7T7tK, KKK decays observed at the dedicated B- meson facilities BaBar and Belle ^7] . In 
particular, in Ref. an experimental evidence for a light and broad scalar resonance 
in the m n7T spectrum of the D — ► -k-r-k decay was found, to be identified with the fo (600), 
and also a less pronounced /o(980) peak. Both the BaBar and Belle collaborations ^7J 
report a distinctive peak around 1 GeV in the m n7T spectrum of three-body 5-decays. 

Our interest in D-decays with an /o(980) in the final state is motivated by two rel- 
ativists quark model calculations of B — > /o(980) and D — > /o(980) transition form 
factors we have performed recently, one of which makes use of the dispersion relation 
approach, the other is being done within Covariant Light-Front Dynamics (CLFD) PQ. 
Both methods are explicitly covariant and require two size parameters that specify the 
unknown /o(980) wave function. In order to deduce these parameters from experiment, 
we fit D(D S ) — > /o(980)X branching fractions, where X can be a pion or a kaon. The 
decay amplitudes are built within the naive QCD factorization approach which neglects 
hard-scattering effects between the spectator and any of the other final quarks. Anni- 
hilation topologies are not considered here. Since m c is less than half the value of m&, 
non-perturbative contributions of order Aqcd/^c may be more important than in B- 
decay amplitudes and the factorization approach thus less reliable. Nonetheless, already 
at tree level, non-leptonic two-body decay channels of D-mesons can be reasonably re- 
produced within this simple factorization hypothesis. This is not the case for two-body 
.B-decays, where penguin contributions are more important. Hence, the decay amplitudes 
are products of a short-distance part, calculated in perturbation theory, and hadronic ma- 
trix elements describing the long distance physics such as the pion or kaon decay constant 
and the D — > /o(980) form factors we are interested in. 

2. DECAY AMPLITUDES IN THE FACTORIZATION APPROACH 

At tree-level order the generic decay amplitudes for a D-meson decaying into a scalar 
/o(980) and a pseudoscalar meson P are given by 



Here, Gp is the Fermi coupling constant, fp is the pseudoscalar decay constant, V cq and 
V uq are the CKM matrix elements with q = d,s and mp> and rrif are the D- and /o(980)- 
meson masses, respectively. Our main interest here is the D — > / (980) transition form 
factor F D ^°(mp). The short distance physics, beyond the scale /i = m c , is described 
by the effective coefficients a\ — C\ + C^/iVc and ct2 = Ci + C%/N c where Ci(fi) and 
C 2 (/i) are the tree-operator Wilson coefficients taken at /i = m c [TH] and N c = 3 is the 
number of colors. The amplitudes for penguin diagrams have the same structure, yet their 
contribution in D-decays is considerably smaller and usually neglected. 

Using the above factorized amplitudes, the decay rates can be computed as 



A(D -> f P) 



Gp 

71 



ajpiml - m 2 fo )F Q D ^ f0 (mp) 
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Table 1 

Experimental branching fractions for D — ► /o(980)X and D s — > /o(980)X (X = 7r,iT) deduced 
from the corresponding three-body decays [8-16] with £>(/o — ► 7r+7r~) = 0.80 ±0.14 and B(fo — > 
K+K-) = 0.11 ± 0.02. JF(D -» / vr ) is the fit fraction of the L>° -> / 7r° decay. 



E791 


B(D + -> / 7T + ) x B(/ -> tt+tt-) = (1.9 ± 0.5) x 10" 4 


FOCUS 
FOCUS 


B(D+ -> / ^ + ) x B(/o -> K+if") = (3.84 ± 0.92) x 10~ 5 
B(D + -> /„#+) x B(/ -> tt+tt-) = (6.12 ± 3.65) x 10~ 5 


ARGUS 

CLEO 

BABAR 


B(D° -> / /?°) x B(/ -> tt+tt-) = (3.2 ± 0.9) x 10~ 3 
£(D° - f K°) x B(/ - tt+tt-) = (2.5+°j) x 10- 3 
B(D° -> /qK ) x B(/ -> #+#") = (1.2 ± 0.9) x 10" 3 


E687 
E791 
FOCUS 
FOCUS 


£(£>+ -> / tt+) x B(/ -> fC+A-) = (4.9 ± 2.3) x 10~ 3 
B(D+ -> / tt+) x B(/ -> tt+tt-) = (5.7 ± 1.7) x 10- 3 
B(D+ -> / tt+) x fl(/ -> tt+tt-) = (9.5 ± 2.7) x 10~ 3 
B(D+ -> / 7T+) x B(/ -> #+#-) = (7.0 ± 1.9) x 10~ 3 


FOCUS 


B(D+ -> /„#+) x fi(/ -> #+#-) = (2.8 ± 1.3) x 10" 4 


CLEO 


i3(D° -> / TT°) ~ B(D° -> TT + TT-TT ) X ^(D° -> / TT°) = 

= (1.1 ± 0.4) x 10~ 2 x (1.0 ± 0.8) x 10- 4 



where |p| = \/[ m z) — (mp + mf ) 2 } [m 2 -, — (mp — mf a ) 2 ]/2m£> is the cm. momentum of 
the decay particles. The transition form factor F Q D ^^°(m 2 p ) is related to the + and — 
components of the pseudoscalar (P) to scalar (S) weak-hadronic matrix element defined 
as (pi is the initial and p 2 the final momentum and q — pi — P2) 

(S(p 2 )\J»\P( Pl )} = ( Pl +p 2 YF + ( q 2 ) + ( Pl - nYF-tf) , J» = qY(l - 7 5 )c , (3) 
by the expression 

In our case Mi and M 2 are the masses of the D- and /o(980)-mesons, respectively. As 
mentioned, computation of the F + (q 2 ) and F_(q 2 ) requires the knowledge of the /o(980) 
wavefunction. In Ref. [I], the /o(980) was taken to be a qq state with a non-strange and a 
strange component which implies a mixing angle. Higher Fock states were not included. 
The /o(980) wavefunction then reads 



fo 



= —f=(\uu) + \dd)j sin6> m + \ss) cos9 m = $ n sin0 m + $ s cos9 m , (5) 

where 8 m is the mixing angle and n = u(d) or s. The spatial component $ n and $ s of the 
wavefunction are of Gaussian form $ n ' s cx exp(— Aa n ^ s k 2 /n 2 ), where k is the modulus of 
the center-of-mass quark momentum and /1 = m q mq/ (m q + m,) = m n>s /2 is the reduced 
mass of the qq pair. The size parameters a nyS are then fitted via the transition form factor 
Fq~* in the decay rates Eq. (J2J) to the experimental branching ratios listed in Table H 
The mixing angle 9 m has been estimated from — ► / (980)tt+ and — > </>tt+ decays to 
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cover the rather wide range 20° < 9 < 40° and 140° < 9 < 160°. There are two solutions 
for 9 m since the mixing angle enters quadratically into the decay rate formula. We choose 
to determine the mixing angle ourselves and thus include 9 m along with a n and a s as a 
parameter, which we determine by fitting the experimental values in Table H 

3. CONCLUSION 

In our effort to model covariant transition form factors for B — > / (980) and D — > 
/o(980) pQ, we have used non-leptonic three-body D decays to fix the parameters a n , a s 
and 9 m of the scalar meson wavefunction. We emphasize that our results are strongly 
dependent on the B(f — > 7r + 7r~) and B(f — > K + K~) fractions and that the uncertain- 
ties of the experimental branching ratios are large (see Table P). Therefore, due to the 
lack of precision and the need for more data, these uncertainties yield large parametric 
correlations and sizable parameter errors. We hope that newer analyses of three-body 
.D-decays with higher statistics will remedy this problem. Definite results for a n>s and 9 m 
as well as B — > / (980) and D — > / (980) transition form factor predictions for space- and 
time-like momenta will be communicated promptly. 
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